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METHOD AND SYSTEM TO PERFORM ENERGY-EXTRACTION 
BASED ACTIVE NOISE CONTROL 

STATEMENT REGARDING 
FEDERALLY SPONSORED RESEARCH AND DEVELOPMENT 

[0001] This invention was made in part with Government support under Grant Number NCC- 
1-01039 awarded by NASA and Grant Numbers 0196198 and 0301740 awarded by the NSF. The 
Government may have certain rights in this invention. 

FIELD OF THE INVENTION 
[0002] This invention pertains to noise suppression, and more particularly, relates to active 
noise control. 

BACKGROUND OF THE INVENTION 

[0003] Sound is a longitudinal wave phenomenon. Interference between waves has been 
known and observed since the time of Sir Isaac Newton. Active suppression of noise has been a 
goal for many decades. The potential applications of active noise suppression are many, ranging 
from aircraft cabins to automobile interiors to medical equipment (e.g., dentist drills, operating 
rooms, etc.) to buildings to household machinery such as vacuum cleaners, lawn mowers, snow 
blowers, washers, dryers, etc. The concept of wave dynamics and the ability of electronic 
amplifiers and moving coil speakers to generate accurate sound waves led to techniques used in 
the 1930s to reduce noise for a harmonic acoustic disturbance. The next important step was taken 
in the 1 950s with the experimental demonstration of a very high-gain static feedback acoustic 
noise reduction scheme. In today's terminology, it is called a collocated control scheme with its 
inherent robustness properties. This basic collocated high-gain static feedback scheme is now 
commercially used for industrial ear protection devices. 

[0004] In the 1970s, the development of acoustic noise control took a slightly different turn. 
Papers such as those by Jessel &Mangiante, and Swinbanks addressed the problem of acoustic 
noise control in ducts. The target application discussed in these papers was the reduction of noise 
in HVAC systems. The main idea was to measure the noise disturbance using a microphone 
located at an upstream location and then to feedforward this disturbance via a controller to a 
speaker located somewhere at the downstream location. The research showed that a proper 
adjustment of the controller gains would arrest the propagation of the disturbance downstream 
from the controlling speaker. This control technique, referred to as the feedforward control, led 
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to several years of active noise cancellation research and development using the feedforward 
approach. 

[0005] For several years, feedforward methods were very popular. It was argued that since 
feedforward scheme doesn't alter system dynamics, there was no risk of destabilizing the system 
(or making it any worse). These schemes required (or otherwise assumed) that the disturbance 
measuring microphone (located upstream) does not sense the output from the control speaker 
(located downstream). If it did, then it would be a feedback scenario with a danger of closed-loop 
instabilities. This meant that the schemes could be successfully implemented only with 
unidirectional microphones and speakers. In general, these conditions are very stringent. 
Nevertheless, researchers came up with ingenious arrangements and obtained impressive results. 
Most well-known are the experiments based on feedforward control by Ross and Roure wherein 
the controller was designed to perfectly cancel the acoustic noise. These experiments involve 
exact knowledge of several system transfer functions, i.e., the design needs to calculate 
differences of transfer functions and invert them leading to high sensitivity of the controllers to 
uncertainty in the system model. Owing to high sensitivity, the model information has to be 
obtained experimentally. This remedy, however, did not alleviate the problems as the controllers 
were very sensitive to the implementation and finite digit arithmetics and had to be implemented 
with on-line adaptation schemes. In spite of these hurdles, their success was a big inspiration for 
much of the work in active control of acoustic noise. 

[0006] However, one important drawback of the prior feedforward strategies that are used is 
the necessity of an on-line adaptation, not due to the nature of the system (acoustic duct) but due 
to the nature of the controller itself. Adaptive control comes at a price. For large bandwidth 
systems such as a 3-D acoustic duct, support of a fast digital signal processor with significant 
memory is essential for the feedforward techniques to work. Therefore, the implementation of 
feedforward controllers really got boosted only after the arrival of cheap, high performance DSP 
(digital signal processor) chips. Another difficulty with feedforward schemes is to ensure that the 
adaptive control algorithm will always converge. Most algorithms guarantee convergence with 
known model, known controller order, and known persistency of excitation but, unfortunately, 
such ideal conditions are very hard to meet. As a result, feedforward methods for noise control 
are generally effective only at specific locations, usually downstream of the disturbance, and at 
specific disturbance frequencies. The need to use high performance DSPs for the 
computationally intensive on-line adaption schemes results in an expensive solution and the lack 
of robustness can make the noise worse in the event of small deviations. 
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[0007] Since the late nineteen eighties, active noise control research took a multifaceted 
approach. During this time, research efforts extended to a multitude of problems, including such 
items as investigation of structure-borne noise, acoustic radiation from structures, reduction of 
noise in ducts and enclosures, and acoustic-structure interaction. Feedforward adaptive 
algorithms work well when a signal strongly correlated with the noise source is available. Most 
of the algorithms converge nicely if the dynamics of the path between the control actuator and the 
sensor is known. Although some earlier successful solutions to structure-borne sound used 
multiple microphones and loudspeakers, subsequent studies showed that loudspeakers alone are 
not very efficient in controlling structure-borne sound. This observation led researchers to invent 
new actuators that directly act on the structure itself. It was found that an effective way to reduce 
the radiated acoustic energy from the structure is to control the structure so as to minimize the 
radiated acoustic power. Although a majority of the research during the 1980's and 1990's still 
used feedforward adaptive algorithms, a large number of researchers were also engaged in 
developing and using feedback techniques. 

[0008] The research in the nineteen nineties especially involved fundamentally new 
approaches to active noise control - a feedback control instead of traditional feedforward 
methods. A significant effort was also expended looking at specific problems (e.g., active noise 
control for aircraft cabins, control of structure-borne noise, and a few others). In spite of 
difficulties, a number of approaches have succeeded in feedback control of reverberant enclosures 
but they are restricted primarily to the low-frequency bands. Additionally, much of the 
experimental work was restricted to small size cavities. Nevertheless, the success stories of 
feedback methods in nineties gave boost to the researchers for further exploration of feedback 
concept for active noise control. 

[0009] The major hurdles in successful implementation of feedforward methods included 
availability of a reference signal having strong correlation to the noise source, necessity of online 
adaptation, and lack of controller robustness. For example, for structure-borne sound in aircraft it 
is difficult to obtain a reference signal which is strongly correlated with the noise producing 
mechanism. This makes feedforward control very ineffective. The advances in robust control 
methodologies coupled with knowledge gained in feedforward control schemes made it possible 
to develop closed-loop stable feedback control methods. Moreover, feedback methods have 
natural robustness compared to feedforward methods which is inherent in the architecture. A 
comparison of LQG (linear-quadratic-Gaussian), rate feedback, and filtered-x LMS (least mean 
squares) algorithm, presented for SISO (single output, single input) feedback control of structure- 
acoustic dynamics, showed that there is a definite rationale to work with feedback methods. 
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[0010] However, designing controllers for acoustic systems is very challenging since there is 
no natural roll-off at high frequencies and the systems are modally very rich. In the mid to late 
nineteen nineties, noteworthy research demonstrating effectiveness of feedback control 
methodologies emerged. For example, a PDE-based (partial differential equation-based) 
feedback control design was employed for controlling acoustic noise in 3-D enclosure using 
piezo-actuated structures. The results were only numerical and the noise was only tonal. 

[0011] Although the feedback control methodology has great potential, it has several 
obstacles to overcome. One problem with many feedback control techniques used is that the 
noise attenuation can be guaranteed only locally (i.e., only at the sensor locations) and not 
uniformly in space and/or frequency. In fact, the noise levels can go high at some other locations 
although there is attenuation at the sensor locations. This essentially is the result of redistribution 
of acoustic energy in space. It is also possible to have worsening of noise levels at the same 
location due to redistribution of acoustic energy in frequency, commonly referred to as the 
waterbed effect. 

BRIEF SUMMARY OF THE INVENTION 

[0012] The invention provides a method and apparatus to provide active noise control to 
reduce noise and vibration in reverberant acoustic enclosures such as aircraft, vehicles, household 
appliances, high precision instruments requiring vibration/noise isolation, industrial equipment 
and the like. The invention allows the designs of quieter aircraft fuselages, vehicles, vacuum 
cleaners, washers, dryers, lawn and garden equipment (e.g., power blowers, lawn mowers, snow 
blowers, etc.), and other industrial/household machinery such as grinding machines, milling 
machines, cooling systems, and the like. The apparatus uses compensators and feedback control 
to render the system "passive" and reduce the structural and acoustic energy in the enclosure. 

[0013] The method to render the apparatus includes the steps described below. A continuous- 
time multi-input multi-output (MIMO) state space mathematical model of the plant (i.e., the 
acoustic enclosure) is obtained. The model is derived via analytical modeling and system 
identification. Compensation is designed to render the mathematical model passive in the sense 
of mathematical system theory. The compensation may be series, feed-forward, feedback, hybrid, 
etc. 

[0014] The compensated system is checked to ensure robustness of the passivity property of 
the plant. The check ensures that the passivity is preserved if the mathematical model parameters 
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are perturbed from nominal values. A passivity-based controller is designed and verified using 
numerical simulations and then tested. 

[0015] Other aspects and advantages of the invention will become more apparent from the 
following detailed description when taken in conjunction with the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0016] FIG. 1 is a block diagram generally illustrating an acoustic environment in which the 
present invention operates; 

[0017] FIG. 2 is a flow chart illustrating the process of creating a controller design in 
accordance with the teachings of the present invention; 

[0018] FIG. 3 is a block diagram of an acoustic enclosure on which the controller design of 
the present invention operates; 

[0019] FIGS. 4a-4e illustrate passification methods according to the teachings of the present 
invention; 

[0020] FIG. 5 illustrates a passivity based controller structure in accordance with the present 
invention for a one-dimensional enclosure system; 

[0021] FIG. 6 is a graph illustrating experimentally obtained frequency-domain noise 
suppression response for the enclosure of FIG. 5 for a broadband disturbance; 
[0022] FIG. 7 is a graph illustrating noise suppression using the teachings of the present 
invention in a three-dimensional reverberant structure for a three-tone disturbance noise at 
frequencies of 100 Hz, 2256 Hz, and 290 Hz; 

[0023] FIG. 8 is a graph illustrating experimentally obtained noise suppression response for 
the system of FIG. 7; 

[0024] FIG. 9 is an isometric view of a three dimensional enclosure having a vibrating 
surface; and 

[0025] FIG. 10 is a graph illustrating estimated and measured frequency responses of a 
controller for the enclosure of FIG. 9. 

DETAILED DESCRIPTION OF THE INVENTION 

[0026] The invention extracts acoustic energy and structural energy in acoustic enclosures. 
The system utilizes acoustic sensors and actuators (e.g., microphones and loudspeakers) and 
actuators (e.g., accelerometers and piezoelectric devices) to reduce the acoustic energy and 
structural energy of the system. A compensator is designed to render the open-loop system 
passive in the sense of mathematical system theory. A positive-real feedback controller is then 
designed. The resulting system is typically a multi-input multi-output feedback control system 
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that incorporates the passifier and the positive-real controller. The system can achieve uniform 
(in space and frequency) broadband reduction in noise and vibration levels in the enclosures. As 
a result, quieter aircraft fuselages, vehicles, vacuum cleaners, washers, dryers, lawn and garden 
equipment (e.g., power blowers, lawn mowers, snow blowers, etc.), and other 
industrial/household machinery such as grinding machines, milling machines, cooling systems, 
and the like can be designed. 

[0027] Turning to the drawings, wherein like reference numerals refer to like elements, the 
invention is illustrated as being implemented in a suitable environment. FIG. 1 illustrates the 
control system 100 and acoustic enclosure 102. The acoustic enclosure 102 may be any of the 
acoustic enclosures described above (e.g., enclosures such as aircraft, vehicles, household 
appliances, high precision instruments requiring vibration/noise isolation, industrial equipment 
and the like. The control system 100 includes acoustic pressure sensors 104 (e.g., microphones) 
and structural vibration sensors 106 (e.g., accelerometers. The acoustic pressure sensors 104 are 
placed within the acoustic enclosure 104 and the structural vibration sensors 106 are placed on the 
walls of the acoustic structure 104. Any number of acoustic pressure sensors and structural 
vibration sensors may be used. The location of the sensors 104, 106 may be application 
dependent. 

[0028] The sensors 104, 106 are in communication with signal conditioners 108, which 
consist primarily of analog filters. The conditioned outputs of the signal conditioners 108 are 
converted to discrete-time digitized signals by analog-to-digital (A/D) converters 110 located 
within energy-extraction controller 112. The energy-extracting controller 1 12 consists of a 
processing unit that stores data and computes the control (or command) signals. In one 
embodiment, the processing unit is a digital signal processor. The energy-extracting controller 
also includes a variety of computer readable media. Computer readable media can be any 
available media that can be accessed by controller 112 and includes both volatile and nonvolatile 
media, removable and non-removable media. By way of example, and not limitation, computer 
readable media may comprise computer storage media and communication media. Computer 
storage media includes both volatile and nonvolatile, removable and non-removable media 
implemented in any method or technology for storage of information such as computer readable 
instructions, data structures, program modules or other data. Computer storage media includes, 
but is not limited to, RAM, ROM, EEPROM, flash memory or other memory technology, CD- 
ROM, digital versatile disks (DVD) or other optical disk storage, magnetic cassettes, magnetic 
tape, magnetic disk storage or other magnetic storage devices, or any other medium which can be 
used to store the desired information and which can be accessed by controller 112. 
Communication media typically embodies computer readable instructions, data structures, 
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program modules or other data in a modulated data signal such as a carrier wave or other 
transport mechanism and includes any information delivery media. By way of example, and not 
limitation, communication media includes wired media such as a wired network or direct-wired 
connection, and wireless media such as acoustic, RF, infrared and other wireless media. 
Combinations of the any of the above should also be included within the scope of computer 
readable media. 

[0029] Digital-to-analog (D/A) converters 1 14 are embedded in the controller 1 12 and 
convert the digital control commands to analog commands. The analog commands are amplified 
by amplifiers 116. The amplified commands are sent to actuators 1 18 to extract energy from the 
acoustic enclosure 102. The actuators 118 consist of one or more acoustic signal generators 120 
(e.g., loudspeakers) at various locations within the acoustic enclosure 102 and/or one or more 
piezo actuators 122 attached to the inside or outside of the enclosure wall. While FIG. 1 shows 
the amplifiers connected to only one acoustic signal generator 120 and one piezo actuator 122, it 
is recognized that the amplifier is in communication with all active acoustic signal generators 120 
and piezo actuators 122. The actuators 118 generate the commanded structural and/or acoustic 
wave motion. Performance monitoring sensors (not shown) may be part of the system to measure 
noise (i.e., pressure) and vibration (e.g., acceleration) at desired locations to monitor the 
performance of the control system. The performance sensors consist of one or more acoustic 
pressure sensors (e.g., microphones) placed within the acoustic enclosure, and/or one or more 
structural vibration sensors (e.g, accelerometers) placed on the walls of the structure. 

[0030] Turning now to FIG. 2, the controller design extracts the acoustic energy and the 
structural energy in the system. The controller design is based on passivity theory. The 
controller is designed according to the following steps. A continuous-time multi-input multi- 
output (MIMO) state-space mathematical model of the acoustic enclosure (henceforth also called 
the "plant") is obtained (step 200). The model is obtained via analytical modeling and system 
identification. Compensation to render the mathematical model "passive" in the sense of 
mathematical system theory is designed (step 202). The compensation may be series, feed- 
forward, feedback, hybrid, etc. (as described herein below). 

[0031] The passivity of the compensated system is checked (step 204) using experimentally 
obtained frequency response data. Numerical tests are performed to check robustness of the 
passivity property of the plant, i.e., to ensure that the passivity is preserved if the mathematical 
model parameters are perturbed from their nominal values. If this is not the case, compensation is 
redesigned (step 206). Steps 204-206 are repeated until the passivity is preserved. A passivity- 
based controller is designed using techniques described herein (step 208). Numerical simulations 
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are performed of the closed-loop system in the presence of a simulated broadband disturbance 
input (step 210). If the performance is not satisfactory, the controller is redesigned and steps 208- 
210 are repeated until performance is satisfactory. 



[0032] Turning back to step 200, it should be noted that many members of industry who are 
engaged in active noise control research resort to experimental identification of system models 
based on frequency response data. A typical procedure involves first obtaining frequency 
response data and then using this data to obtain necessary transfer function or state-space models. 
These identification techniques also have limitations as the currently available system 
identification algorithms work well only for low-order systems and systems with nicely separated 
modes. For a 3-D acoustic enclosure with a vibrating boundary surface, these identification 
techniques face difficulty even with two-input one-output system. The invention provides a 
general purpose analytical modeling methodology for a 3-D acoustic enclosure having a vibrating 
boundary surface. This modeling technique accounts for the acoustic-structure interaction 
dynamics and yields time-domain as well as frequency-domain models. 

[0033] A non-homogeneous wave equation is derived below for a general 3-D enclosure 
wherein the source term in the equation represents the acoustic disturbance arising from the 
vibrating plate boundary. Since a closed-form analytical solution is not possible when the entire 
boundary is considered to be vibrating, an assumption is made where only part of the boundary is 
considered to be vibrating. First, an infinite dimensional model is obtained which accounts for 
the structural-acoustic coupling. Next, the assumed modes method is employed to obtain a finite- 
dimensional approximate model, which is used for control design purposes. Both the state-space 
model and the transfer function model are derived analytically. 



[0034] Turning now to FIG. 3, a 3-D acoustic enclosure 300 with a part 302 of one boundary 
surface 304 representing a vibrating plate (denoted by S\). Let the entire boundary surface be 
denoted by S. Prior to proceeding, given below is the list of definitions of various variables and 
parameters used in deriving the acoustic-structure interaction dynamics for the enclosure. 



v(/, x, y, z) is fluid velocity in m/s 

pit, x, v, z) is acoustic pressure fluctuation from equilibrium in N/m 2 

p is fluid density fluctuation from equilibrium in kg/m 3 

po is acoustic pressure at equilibrium in N/m 2 

po is fluid density at equilibrium, kg/m 3 

© (/, x l , v, ) is the transverse displacement of the plate in m 

f (J, x, , v, ) is the externally applied force on the plate due to acoustic pressure from an 
exterior noise field in N 
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«,<0 are control forces acting at the points (*/,; >>/,); /" = 1 , 2, . . . , r on the plate and 

normal to the plate in N 

[0035] The acoustic pressure and the fluid density at any instant of time are given by p(t, x, y, 
z) + p 0 and p(t, x, y z) + p 0 . The control forces, «,(/), are produced by piezoelectric actuators and 
can be treated as point sources. For theoretical derivation, we make a necessary and reasonable 
assumption that the plate is simply supported on its four sides and that the initial displacements 
and velocities at all points on the plate are zero. The goal is to obtain a governing relation 
between the acoustic pressure fluctuation p(t, x, y z) inside of a 3-D enclosure such as enclosure 
300 in FIG. 3 and external disturbances / (/, x x , y, ) and uff), i = 1 , 2, . . . , r. The model of the 
vibrating boundary surface (S\) is first established followed by the development of the acoustic 
model for inside of the enclosure. 



[0036] The dynamic model of the vibrating surface 304 (i.e., the plate) can be obtained using 
Hamilton's principle, which states that: Of all possible paths of motion that can be taken between 
the instants of time t\ and t 2 , the actual path taken by the system gives a stationary (extremum) 
value to the integral 

I=\\K X +W x +U x )dt (1) 

where, K\ is the kinetic energy of the system, U\ the potential energy of the system, and W\ is the 
work of the externally applied forces. It is assumed that the quantity K\+Wi- U\ is conservative. 
Now, for a two dimensional system with one dependent variable, if the integral in Equation 1 can 

be expressed as 

h 

K<o) = j$jF(x x ,y x ,t,cD,a>,co xX ,G) yX ,G) XiXi y co Xyyx )dx x dy x dt (2) 
/, s, 

from Hamilton's Principle and variational calculus, the governing PDE (partial differential 
equation) for the system can be derived from 



where 



e x (a>)-e 2 (<6) = 0 (3) 

dF w dFa d'F^ d'F^ d'Fa 

£ l (a>) = F ^ — + ™ 

o> dx x dy x dx x 2 dy\ dx x dy x 

dF. 
£ 2 (d>) = - 



dt 

[0037] For a vibrating plate, one can calculate U\ 9 K\ 9 and W\ as follows: 
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1 2 J ~ y P K dt ' 1 
6 1 



(4) 



'-' s i 

where /», and p p are the thickness and density of the plate, respectively, E is the Young's modulus, 
v is the Poisson's ratio, and is the flexural rigidity of the plate. In Equation 2, the 

12(l^ 2 ) 

function F is given by 



F = -£{(S + !f)-2(l-v)[. 



d 2 wd 2 co d 2 m 2 h ,3<» 
2"dx 2 ' dy 2 ' " v * '' l dx 2 dy\ K dx x y, } U+ 2 Pp( ~di' ) 



+ £ u, (t)a(t, x, , y, )S(x, - x u )S(y, - y u ) + [f(t , *,,>>,) + x, , jy, ,0)]to (r , x, , .y, ) 



(5) 



[0038] The governing equation of the plate can be obtained from Equation 3 as 

Z>V 4 » + = -x^C, -* ( ) + /(',*,,J'i) + P(',*i>.M) (6) 



[0039] The boundary conditions for the plate are given by 



d'a>(t,0,y) _ 



= 0 



^(f,0,v) = 0, 

<o(t, jc,0) = 0, 
= 0, 



= 0 



d z a>(t, X ,Q) _ q 
d 2 a>(t,xj y ) 

* — - — u 



(7) 



[0040] The governing equation for the acoustics in the enclosure 400 shall now be derived. 
For the enclosure 300 under consideration, a part of one of its boundary surfaces is a vibrating 
surface and we can consider this surface as a distributed noise source. The fluctuating volume 
flow per unit volume is then given by 
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p 0 co t (t, x, , y x )5 (z) for 0<x\<b-a and 0<yi<d-c 

[0041] The linearized mass conservation equation and momentum conservation equation for 
the system under consideration can be written as 

dp(t,x,y,z) + ^ v ^ ^ z) = pQ0}( ft ^ ^ w (g) 

A> a ^y >Z) + Vp(r, W ) = 0 (9) 

[0042] Now, using the relationship p(t, x, y, z) = c 0 2 p (t, x, y, z) , where c 0 is the sound speed, 
the divergence ( V •) of the momentum conservation equation, and then subtracting the time 
derivative of the mass conservation equation, the non-homogeneous wave equation is obtained as 
follows 

v 2 P (t,x,y,z)-\ d2p(t : x ; y > 2 K - Po ^ 

c 0 dt (10) 

H(x - b)][H(y - c) - H(y - rf)]S (z) 

where H(») is the Heaviside function. For the closed 3-D enclosure 300 under consideration we 
have the following boundary conditions 

v(t,x 9 y,z)-n = 0 on S (11) 

where n is the outward normal direction at every point on surface S. Using Equation 9, the 
boundary conditions in Eq. (11) can be rewritten as 

Vp(t,x,y,z)-n = 0 on S (12) 

[0043] We will now derive the governing differential equations for the vibrating surface 
(plate) as well as the acoustic pressure fluctuation dynamics in the modal space as it will allow us 
to obtain the finite-dimensional approximate models of the system for controller design and 
experimental verification. First, we will consider only the homogeneous part of Equation 6 with 
boundary conditions (i.e., Equation 7), that is, 

DV 4 co + hp p ^- = 0 (13) 

[0044] Assume the solution is of the form co(t 9 x l , >>, ) = Rc{W(x ] , y x )e it0 ' } . Equation 1 3 is 
transformed into 
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V A W -(K&yw = 0 (14) 
where K? = hp p/ D, Q 4 =co 2 . Now separating variables x\ and y\, W can be written as 

W(jc,j; 1 ) = X(x 1 )Y(y,) 
Equation 14 can now be re-written as 

[X (4) -(KnyX]Y + 2X n Y" +XY i4) = 0 (15) 
[0045] Now considering the solutions which are separable. Let 

X" = -$ 2 X (16) 

then 

X i4) = -p 2 x" = p 4 x 

Similarly, the equation to determine Y is 

r (4) -2p 2 Y" -[(KG)* - p 4 ]7 = 0 (17) 

[0046] For the case, (ATQ) 4 - p 4 > 0 , the solutions of Equations 16 and 17 are given by 



X(x x ) = A x cos(/fc 1 ) + A 2 sin(pjt! ) 

Y( v, ) = B x cos(^, ) + B 2 s'm(Sy x ) + B 3 cosh(^ ) + B 4 sinn^ ) 



(18) 
(19) 



where S = ^j(KCi) 2 - p 2 , s = ^(KO) 2 + p 2 , and A iy 5/ are coefficients to be determined. Note 
that for the case when (ATQ) 4 - p 4 < 0 , there are no solutions. The solution W{x\j\) can now be 
given by 



W(x x , y l ) = cos( Ai ) + ^ 2 sin( J3x x )] x 

[5, cos(<9y, ) + B 2 sin(<5y, ) + B 3 cosh(£y, ) + B 4 sinh(sy, )] 



(20) 



[0047] The boundary conditions (7) can be expressed in terms of W as 



W(0,y x ) = 0, 
W(b-a,y x ) = 0, 
W(x l ,0) = 0, 



d'iV(0,y x ) _ 



= 0 



d 2 fV(b-a,y x ) _Q 



9V(^,0) _q 
d 2 W(x x ,d-c) 



(21) 



= 0 
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[0048] The coefficients and parameters in Equation 20 can be determined from these 
boundary conditions to be A\ = B\ = B 3 = B 4 = 0, sin[(6 - a)P] = 0, sin[(</- c)8] = 0. The normal 
modes of the plate are 

<phm(xi, Vi) = Anm sin(y&xi) sin(5«y, ), n, m = 1 ,2,— (22) 

where, pn = nn /(b -a),Sm- mn l(d - c). The natural frequencies are 

«u =«L =^[(ir-) 2 + (^-) 2 ] (23) 

A o — a a — c 
By the orthogonality condition: 

the coefficients can be determined as A nm — 2. 

[0049] The excitation of the plate due to acoustic pressure load from inside of enclosure is 
insignificant compared to piezo forces. That is, the dynamic interaction is much stronger from 
structure to acoustics than from acoustics to structure. Data obtained in applications such as 
aircraft cabins shows that the dynamic interaction from fuselage to interior acoustics is much 
stronger than in the opposite direction. In view of this, the acoustic pressure fluctuation p(t, x, y t 
z) term in equation 6 is neglected without much loss of accuracy. The mathematical simplicity 
achieved by this assumption outweighs the effects of modeling errors caused by the omission of 
this term. As a result, the governing equation 6 for the plate can now be re-written as 

DVWA/7,^ = 5>i(0*(*i -x Xi *iy x -y v ) + f{t,x uyi ) (24) 
eft ,=i 

[0050] Using the assumed modes approach the solution of this PDE can be expressed as an 
infinite series representing an infinite sum of the weighted normal modes of the system, L e., 

OO O0 

<o(t, x x , y x ) = x Z (0*~ Oi > y\ ) ( 25 ) 

n=\ m=\ 

where q nm (f) are the modal coordinates and 4> IBW (x 1 , i y 1 ) are the corresponding mode shapes. For 
the control design purposes, we need to obtain a finite dimensional model of the system. This can 
be done by approximating the solution in equation 25 using a finite number of modes. That is, we 
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can truncate the series by choosing limits on n and m as n : 1 — > p x and m : 1 -> p 2 . Then, 
the solution for co (•) can be given by 

<o(t,x x y y x ) m £ ^q nm (tys> nm (x l9 y, ) (26) 

»=1 m=l 

[0051] Substituting equation 26 into equation 24 and multiplying it by ty nm (jCj , y x ) and 
integrating it over Si, we get 

hp p 4„ m {t) + D(P 2 n +S 2 m ) 2 q m (t) = ±-£u,(M m (x u ,y ll ) + f m ,(t) (27) 

where, 

f nm (') = -5- J J/0 , x, , * » _ (x, , V, )</S, 

Si s , 

[0052] It should be noted that equation 27 does not contain any damping terms as there is no 
known method to analytically derive the modal damping in the system. The damping has to be 
introduced a posteriori into the model. Introducing the modal damping 
t^„ m (n = 1,2,..., p l9 m- 1,2,..., p 2 ) into the equation, we obtain 

q nm (0 + IC^J^ (/) + co 2 nm q nm (0 = J- [-1 £ w . (04> Bflf (*, , v, ) + /_ (/)] (28) 

The modal coordinates can be determined by Duhamal integral. Once the modal coordinate 
solution is obtained, the response of the plate to external and control forces can be obtained by 
equation 26. 

[0053] The solution for acoustic pressure /?(•) using the assumed modes method. Let us 
assume that the solution of the wave equation 1 0 with boundary conditions 1 2 can be expressed 
as a finite sum of weighted normal modes. 

1 h h 

P(t,x,y,z) = X Z IXm3 ( *'^' Z) (29) 

*, =0A 2 =0*3=0 

where V*,*^ (x,y,z) are the normal modes of the homogeneous part of equation 10 with 
boundary conditions (equation 12). That is, 

¥V 2 * 3 (x,y,z) = A kxHh cos^xcos^cos^z (*i,fe,*3 = 0,1,2,...) (30) 
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where ^ = k x n ll x , ^ = * 2 7t ll y , ^ = k 3 n ll z and q^k^ (0 316 me modal coordinates. 
Substituting equations 26 and 29 into wave equation 10, we get 




n=\ m=\ 



Multiplying equation 31 by V*,*^ (x 9 y,z), integrating the equation over entire volume of the 

enclosure 400, using the orthogonality condition, and adding the modal damping term to the 
equation, we obtain 

where ^ are modal damping ratios, 



[0054] Now that the solution for acoustic pressure has been described, the state-space models 
of the system shall be described. Let the state variables be defined as 

x = {x\x°y = [(*L,yL) T ,(x a kik2ki ,y\ Hh YY (33) 

where x p nm = q tm9 y^ m = q„ m ,x a =Q kkk >y a kkk = <L hhh > n =1,2,.. .pi; m = 1,2,.../^; 

*1*2*3 *1*2*3 *1*2*3 *1*2*3 

^1—0,1,. fe=0,l,.. .,/2; and £ 3 =0,1,.. .,/ 3 . If the external force can be expressed as 
/('> *i > >>i ) = f(f)g(x x ,>>,), f nm (0 can be expressed as 

/„*,(') = [4" J f^iJ,¥ I .(x,^>fi I l/(0 = Y«/(0 (34) 



Then equations 28 and 32 can be re- written as 
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x p = v p 

nm Stun 

Jf" = V° 
*1*2*3 -^*1*2*3 

Aft f 2^ ^ 2 a 

<" ~~ F~ a A,* 2 A 3 « m )>'»« + > , * 1 * 2 A 3 = _fi, *i* 2 *3 X *l*2*3 ~ ^^*1*2*3 C ° *1*2*3 ^*1*2*3 



(35) 



[0055] In the state-space form, the model for a 3-D enclosure under consideration is given by 
Ex(t) = + Bu(t) + Df(t) (36) 



Matrices E, A, B, and D are given by 



E = 



E u 0 
L E 2l E 22 



A = 



0 A 



22 _ 



B = 



hp 0 S, 



0 



hp 0 



D n 
0 



where, En = I and An = diag( A™ ) are square matrices of order p\p 2 , E 22 = I and A22 = 
diag( A*? 2 * 3 ) are square matrices of order (/1 + l)(/ 2 + l)(/ 3 + 1), Bn is a p { p 2 x r matrix, and 
Dn is a p\p 2 x 1 matrix. Moreover, matrices E 2 i, An, A 22 , Bn, and Dn are given as follows 



E 2l = 



c oPo 



0 0 —0 0 

0 ^oom a oo\p x p 2 

0 0 • "0 0 

0 a W3» '" 0 "VA^ 
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0 1 



0 1 

,2 



0 

<^,/>2 Oll'-Vll) 



Al = 



l7 p i p 2 



[0056] If we place k microphones in the position of (x f , y. , z, ), / = 1, . . . , k , the outputs of the 
system can be written as P(t ) = (/?(/, x, , , z, ),..., p(t, x k , y k , z k )) T . Since matrix E is invertible, 
the state-space model (equation 36) can be re-written as 

x(t) = E~ x Ax(t) + E- l Bu(f) + E- l Df(t) 
P(t) = Cx(t) 
where 

c=[c; c 2 ... c A f 

The C matrix has £ rows, and C, has the form 



C,. = (0, ■ • • ,0; if/ oox ( Xi y. , z, ),0, - - • ^ (*.>;,. , z,- ),0) i = 1, • • • , k 



(38) 
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[0057] Now that the state space model has been described, the transfer function model shall 
be described. For single-input single-output (SISO) systems, if we let /(/,*, ,>>,) = 0 and take the 
Laplace transform of equations 28 and 32, we get 

(s 2 +2^^ + fi ii)a„W = r ^<|)^(x n ^„)f/W (39) 

h P P s \ 



2 P\ P2 

(s 2 +^ i¥ « ltt n^ ¥ pa tt W = ^II« W) ^ 2 2-W (40) 



1a * ' 

that is 



Qtfyk. CO = 2 5 ZZ 2 t • ^O) 



(41) 



where J - • Taking the Laplace transform of equation 29, we can obtain a SISO transfer 

function for the output being the acoustic pressure measured at the microphone located at 

(x 0 ,y 0 ,z 0 ) as 

[0058] To obtain the transfer function of the enclosure (i.e., equation 42), we used "hard" 
boundary conditions (i.e., equation 1 1 or equation 12. For more general boundary conditions, 
such as the soft boundaries, the following transfer function is used. On the boundary S, there 
exists the following relationship 

P(r, s) = Z(r, s)U(r, s) - n(r) (43) 

where r = (x,y,z), P(r,s) and U(s) are the Laplace transforms ofp (/ , r) and v(/,r) respectively, 
and Z (r ,s) is the wall impedance. Substituting equation 43 into the Laplace transform of 
equation 9, we can get 

VP(r,j)-*(r) = -p 0 s ffi^ where r in S (44) 

[0059] Assume that the solution of the Laplace transformed equation of the wave equation 
(equation 10) can be written as 
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^■ZWffiW (45) 

i=i 

where the following orthogonality condition are imposed on mode functions 

\%(r)V J ( r )dV=& l) (46) 

V 

[0060] Multiplying the Laplace transform of the wave equation equation 10 by <£ y (r) and 
take the volume integral of the equation to get 

j¥,(r)V 2 P(r, -4 \^j(r)P{r 9 s)dV = -p^f^^Q^s) (47) 

V ^0 V n=\ m=l 

where rj Jlm = J«P J (r^ rall (x-a^-c)[//(x-a)-//(x-6)][//(j,-c)-//(^- d)^(z)dV . Using 

Green's first formula in terms of the appropriate variables, the first term in the left hand side of 
equation .47 can be expressed as 

jV. (r)V 2 />(r, s)dV = jV, (r)VP(r, s) • n(r)dS - jW, (r) • VP(r, s)dV (48) 

r s K 

[0061] We can then calculate 

jT,(r)VP(r,s) -n(r)dS = - Po *Z[^, (*)£(*)] 

s /=i 

JV^ 7 (r) • Vi>(r , = X P„a W 

K /=1 

J^(r)/»(r^) < /K = e,( S ) 

V 

where D„(s) = J T ^^ (5) ^S, p, 7 = J^.(r)V^(r)^. Equation 47 becomes 

S ' V 

X[Po^W + ^JaW + 4aW = A^ 2 i;Zny nra a m W y=U,-,/ (49) 



[0062] Note that in practice it is often found that Dji(s) for j ^ i are negligible, in that case the 
expression for Q/s) can be written as 

QjW~- — — Z&^e„„w (so) 

s - + Po sD Jj {s) + p ji 
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[0063] Substituting equations 50 and 39 into equation 45, we get the transfer function for the 
enclosure with soft boundaries as 

G (j ^ = y PbJ 2 gp %M r f^ yjj^j 1 rsn 



[0064] Now that the modeling has been described, the linearized dynamics of a structure- 
acoustic enclosure can be expressed in the state-space form as: 

x(t) = Ax(t) + Bu(t) + B N N(t) (52) 
y(t) = Cx(t) + Du(t) + D N N(t) (53) 
z(t) = C 2 x(t) (54) 

where *, «, # and z are w- w- and 1- dimensional state vector, input vector, control-output 
vector, and performance-output vector, respectively (all functions of time /). The state vector x(t) 
consists of the acoustic state variables (pressures and rates) as well as structural vibration state 
variables (modal amplitudes and rates for the elastic modes). N(t) is the external disturbance 
vector consisting of external noise and/or forces (e.g., due to the aircraft engine or air-friction at 
the airframe). The control-output y(t) is utilized for generating the feedback control signal u(t), 
while the performance output z(t) denotes a combination of noise levels (pressures) and vibration 
(accelerations) at locations where the performance is monitored. A, B, C, D, B M D N are 
appropriately dimensioned matrices. The relationship between u(t) and y (t) can also be expressed 
as the transfer function 

P(s) = C(sl - Ay 1 B + D (55) 

[0065] The controller design is based on the concept of passivity. An mxm rational matrix 
G(s) is said to be positive real (PR) if 

(i) all elements of G(s) are analytic in Re(s) > 0; 

(ii) G(s) + G*(s) > 0 in Re(s) > 0 (where G* denotes the conjugate-transpose); or 
equivalently, 

(iia) poles on the imaginary axis are simple and have nonnegative-deflnite 

residues, and 

(iib) G(jco) + G\j(o)> 0 for <d g (-00,00) 

For single-input, single-output (SISO) systems, condition (iib) is equivalent to: 
Re[G(_/o>)]> 0 , or equivalently, the phase of G(ja>) remains between -90° and +90°. 
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[0066] Gfs) is said to be weakly strictly positive real (WSPR) if G(s) is stable (i.e., has all the 
Smith-McMillan poles in the open left-half of the complex plane), and the inequality in (iib) is 
strict (>), which requires the transmission zeros to be in the open left-half plane. The stability 
requirement can be weakened as follows to define the least restrictive class of positive real 
matrices, which permits imaginary-axis poles. 

[0067] Aximxm rational matrix G(s) is said to be marginally strictly positive real (MSPR) if 
it is positive real, and 

G{jcoi) + G*(jco) > 0 for G> € (-00,00) 
Note that if G(s) is stable and MSPR, it is WSPR. 

[0068] The most important property of positive real systems is that they are robustly 
stabilized by any MSPR controller C(s) if none of the purely imaginary poles of G(s) is a 
transmission zero of C(s). Thus a PR system is robustly stabilized by any WSPR controller. The 
stability is guaranteed regardless of modeling inaccuracies and uncertainties. 

[0069] This robust stability property offers an attractive approach to robust controller 
synthesis for PR systems that can provide optimal performance along with guaranteed stability. 
However, the transfer function of an acoustic-structure system is not positive-real, (i.e., it is a 
non-passive system). The approach of this invention is to "passify" (i.e., render passive) the 
system using compensation. Subsequently the feedback loop is closed with a suitably designed 
positive-real controller. 

[0070] The passifying compensators include series compensation, output feedback 
compensation, feedforward compensation, hybrid compensation (i.e., combination of series, 
feedback, and feedforward), and sensor blending and/or actuator allocation (applicable to systems 
having multiple sensors and actuators). FIGs. 4a-e shows block diagrams of different types of 
passifications. FIG. 4a illustrates series passification. FIG. 4b illustrates feedback passification. 
FIG. 4c illustrates feedforward passification. FIG. 4d illustrates hybrid passification. FIG. 4e 
illustrates sensor/actuator blending. In each of FIGs. 4a-4e, P(s) indicates the non-passive plant 
and P c (s) indicates the passified plant. FIG. 4-e shows one of many possible configurations for 
hybrid passification. For single-input, single output (SISO) systems, passification can be 
accomplished in a fairly simple manner by using Bode plots, whereas for multi-input multi-output 
(MIMO) systems, passification is more complex and is best performed in state-space setting using 
tools such as linear matrix inequalities. 
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[0071] Different types of non-passive systems warrant different types or combination of 
different types of passification methods. An important consideration in the design of a passifying 
compensator is the robustness of passification. The stability robustness depends on the 
robustness of passification as the passivity is not inherent to the system, but depends on the 
passifying compensator. The structural-acoustic system is usually non-minimum-phase. 
Therefore the approach of this invention is to first render the system minimum-phase by using a 
feedforward compensation Cff(s). For many structure-acoustic systems, a constant feedforward 
matrix Dffis sufficient. The design of such a matrix can be accomplished by using linear-matrix- 
inequalities (LMI). Only the minimum necessary feedforward gain should be used for 
passification because it tends to limit the effective feedback gain and hence the performance. 

[0072] Designing a series passifying compensator can be accomplished by using LMI 
techniques. Alternatively, if large numbers of actuators and control sensors are available, the 
techniques of optimal blending of sensor signals and optimal allocation of control inputs are used. 
After designing the passifying compensator, the passivity of the compensated plant is verified 
using experimentally obtained frequency response data of the system and numerically generated 
frequency response data of the passifying compensator. In addition, numerical computations are 
also performed to verify that the system transfer function remains 'robustly PR 1 in spite of 
parameter variations within a specified range (such as +/- 25%). If the system does not remain 
'robustly PR', the passification is redesigned to accomplish robust positive realness. 

[0073] Designing an optimal controller that is also PR. The methods described in the paper 
"Synthesis of LQ-Optimal Constant-Gain Positive-Real Controllers", Control and Intelligent 
Systems, 29(3):65-73, 2001 by A.G. Kelkar, Y. Mao, and S.M. Joshi, hereby incorporated by 
reference in its entirety, for optimal constant-gain positive real controllers are used. Closed-loop 
simulation studies are performed, using a white noise disturbance input signal co(/), to determine 
the reduction in noise and vibration. If the reduction is not adequate, an optimal PR linear- 
quadratic-Gaussian (LQG) controller is designed. The performance is verified via simulation. If 
the performance is not adequate, additional actuators and control sensors are installed and the 
design process is repeated. If the performance is satisfactory, the overall controller, consisting of 
the passifying compensator and the PR feedback controller, is realized as a digital filter and is 
implemented as the "energy-extracting controller" shown in Fig. 1 . 

[0074] Now that the theory has been described, the details of the passification techniques 
shall be described. Passification can be performed via series compensation. Series compensation 
can be performed either at the input (pre-compensation), as shown in FIG. 4(a), or at the output 
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(post-compensation). We shall only present the pre-compensation technique. The series post- 
compensation is similar. 

[0075] Consider a non-passive plant: 



Series passification can be used only for Lyapunov-stable, minimum-phase systems. Acoustic- 
structure systems are always asymptotically stable, but they can be non-minimum-phase (i.e., one 
or more transmission zeros of the realization [A, B, C, D] can be in the closed right half of the 
complex plain. If the system is non-minimum-phase, it must first be rendered minimum phase via 
feedforward compensation. If an m x m constant feedforward gain matrix Z^is used, the 
resulting system is [A, B, C, (D + Dj$\, and its inverse is given by: 

[(A - BD;jC),BD;j,D;jC,D;j]v/hQTe D eff = D + D ff is assumed to be nonsingular. Therefore the 
transmission zeros of [A,B,C,(P + are the eigenvalues of (A - BD'^C) . Thus, if the system 
[A, B, C, D] can be stabilized by constant-gain output feedback F, it can be rendered minimum- 
phase by a constant feedforward gain gain Dff which is obtained from: F= (D + Djff 1 . F can be 
obtained by using multivariate pole placement methods or by trial and error by setting Dp to be a 
positive diagonal matrix and increasing its elements until the inverse system is stable. It is 
advisable to find the smallest £^such that the inverse system is stable. The resulting minimum- 
phase system is henceforth denoted as: [A,B,C,D]. 

[0076] The series passifier (Fig. 4(a)) is given by: 




(56) 




(57) 



Then the compensated system is given by: 




(58) 



where 




,B = 



(59) 
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For the compensated system (G c (s)) to be passive, we need ( A , B , C , D ) to satisfy the following 
inequality 



A T P + PA PB 
B T P 0 









T 












c 


D 




0 


-/ 


c 


D 


+ 


















0 


I 




-/ 


0 


0 


I 



<0 



(60) 



[0077] Matrix P is a symmetric positive definite solution to the inequality (i.e., Equation 60). 
Existence of such P ensures positive realness of Gc(s). Note that Equation 60 is nonlinear in 
matrix variables Ac, B c , C c , D c , and P, and therefore, does not constitute an LMI. This nonlinear 
matrix inequality, however, can be converted into an equivalent LMI using a judiciously chosen 
nonlinear matrix transformation. Given below is the procedure to convert this nonlinear matrix 
inequality into an LMI. 



[0078] Let us suppose matrices P and P can be partitioned as: 

'X 



TY Nl JX Ml 

[n t * J [m t * J 



where X and Tare symmetric and full rank matrices with the restriction that 

XY + MN 1 ' =1 
Further, if we define 

'X I 



M T 0 
l _ 



n 2 = 



/ Y 

0 N T 



then, since PP" =1, II i and II 2 have the following properties 

pu x =n 2 and nfpn, =nfn 2 = 



X I 
I Y 



(61) 



(62) 



(63) 



(64) 



Let n = I ^ . Now pre- and post-multiplying Equation 60 by Yl and Yl, respectively, does 
not change the inequality, i.e., 



n T x (A T P + PA)Tl, U]PB 



(B T P)Il x 



0 



0 



-n\c 7 



-CTl, ~(D + D r ) 



<0 



(65) 
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where 



7 oJa bc c 1 

" Y n\q A c 



X I 
M T 0 



AX + BC C M T Al 
YAX + YBCM T + NAM T YA\ 



[0079] Now consider the following change of variables 

A = YAX + YBC C M T + NA C M T 
B = YBD C + NB C 
C = C C M T 
D = D C 

Using transformation equations (66), we get 



(66) 



nl2 T pn l = 



AX 



+ BC A 1 
A' YAJ 

\(AX + BC) T A T 

(YA) T 



L ^ (YA)^ 

\y n\[b c J [b 
cn,=[c dc^ X mT lyicx + ncc] 



(67) 



Substituting Equations 67 in Equation 65, we get an LMI in variables X r Y f A,B,C, andZ) 



A** 

A + A T 

D T B T -CX -DC 



(*) (*) 
YA + A T Y (*) 
B T -C £>** 



<0 



(68) 
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where A** = AX + XA 7 + BC + C T B T , D** = -DD- D T D T and (*) indicates the transposes of 
their corresponding off-diagonal counterparts. 

[0080] If the solution to Equation 68 is feasible, the passifying compensator Cs(s) can be 
recovered by following the step-by-step procedure below. 

1 ) Solve equation 68 to obtain X Y, A, B, C, and D 

2) Construct M, N, and P such that they satisfy Equation 64. (Note that M and N should 
satisfy Equation 62. Also, since ^ > 0 , Y>0, and / — XYis nonsingular, i.e., we 

can always find nonsingular square matrices M and N satisfying equation 62) 

3) Solve equations 66 in reverse order to obtain Ac, Be, Cc, and Dc. 

[0081] Passification via feedforward (parallel) compensation shall now be described. 
Consider a plant: 




(69) 



and feedforward compensator (FIG. 4 (c)): 




(70) 



Then the compensated system is given by: 
_ \x = Ax + Bu 

y = Cx + Du 



(71) 



where 




and 



(72) 



C = 




For the compensated system to be passive, we need, as in the previous cases, 
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B T P 











T 










PB 


+ 


C 


D 




0 


-/ 


c 


D 


0 




0 


I 




-/ 


0 


0 


I 



<0 



(73) 



Pre- and post-multiplying equation 72 by n T and IL respectively, we get 



n T x (A T P + PA)Tl x U\PB 
(B T P)Yl l 0 
where 

i\[pau x = n T 2 2n x = 



0 -Y1\C T 
-CTI 1 -(D + D T ) 



AX A ' 

YAX + NAM 7 YA _ 



<0 



(73) 



(74) 



[0082] Now consider the following change of variables 

A = YAX + NA C M T 
B = NB C 
C = C C M T 
D = D C 
Using equations 75, we get 



(75) 



n l T A T pn l = 



AX A 1 



nfp/? = n^ = 



5 



cn, = [cx + c c] 



(76) 



Substituting equations 76 in equation 73 we get an LMI in variables X, Y, A,B, C, and D . 
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(*) (*) 
YA + A T Y (*) 
B T Y + B T -C D L 



<0 



(77) 



where D ± = -(D + D 7 + D + D 7 ) and (*) represents elements which are transpose of their 
corresponding off-diagonal counterparts. The design process will again involve the steps given in 
previous paragraphs. In this case also, similar to series passification. it is required that the non- 
passive plant be open-loop stable. 



[0083] Constant-gain feedforward passification shall now be described. The simplest 
feedforward compensation is a constant gain Dc rather than a dynamic compensator. In that case, 
the compensated system is PR if and only if there exist matrices P=P T > 0 and Dc such that 



A T P + PA 



PB-C T 



B Tp_ c - (D + Dc) _ (D + Dc y 



c) J 



<0 



(78) 



Dc is constrained to be symmetric and non-negative definite and a solution that minimizes the 
trace of Dc is obtained. 



[0084] Passification via sensor blending and control allocation shall now be described. In 
some cases, many sensors (accelerometers and/or microphones) may be available for measuring 
the vibration accelerations and pressures at various locations. In such cases, sensor blending can 
be performed (matrix M in Fig. 4(e)) to passify the system. Similarly, if many actuators (piezos 
and/or loud speakers) are available, control allocation can be performed (matrix N in Fig. 4(e)). 
These passification techniques for obtaining a system that is robustly positive real in the presence 
of parameter uncertainties are described herein. 

[0085] Consider the system S (p) given by: 

x = A(p)x + B(p)u; y = C(p)x + D(p)u (79) 

where pe R k denotes the vector of uncertain parameters and x, u, y are «-,r-,and /-dimensional 
state, input, and output vectors respectively. A(p) y B(p), C(p), D(p) are appropriately 
dimensioned matrices that are assumed to be affine functions of the parameter vector p. The 
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parameters p t (i =1 ,2 , ■ • £),which are components of p, are assumed to lie in a hyper- 
rectangular box in the parameter space, defined as: 

P ={p'Pi^ [£,,¥,], i = 1,2,...*} (80) 
/?, and Pi , / =1 ,2 , —k ,denote the lower and upper bounds on p h Define 

(81) 



= ^(/>i+P,.)and8,. = -p.) , 7=1,2, . . .k 



[0086] Expanding p t as p f = Mi + afii where a, e [-1 ,1], the system matrices can be written 



^(p) = A( M ) + 2>,M/ ; = 5(/i) + l>5,i?, 

1=1 7=1 

C(p) = C( M ) + Jar^C,. ; Z)(/7) = /)(//) + X a > 6 /A 



(82) 
(83) 



where // = [//i, M2,'"Mk] and A/, Bj, Q, D, are constant matrices. The fth element 
(i =1,2, —A) of the yth vertex r/ , (J =1,2, -2* )of £ is given by: 

Pi = Mi (Yy = -lorl) 

(84) 

[0087] Then the system matrix at/? =// is given by 
^(/>0 = ^(//) + I>A4 



B(//), C(//), and D(j/) are similarly defined. 

[0088] The following theorem presents a sufficient condition for the system to be robustly 
passive V p e B. Note that passivity is defined only for 'square 'systems (i.e., with the same 
number of inputs and outputs). Suppose / = r. The system S(p) is passive if there exists a 
symmetric positive-definite matrix P that is a solution of the following system of 2* LMIs: 
~A T (pj)P + PA(p J ) PB(p J )-C T (p J ) 
_B T (p J )P - C{p j ) - (D(p J ) + D T (p J ))J 



Z(p\P) = 



£0 for/ =1,2, -2* 



(85) 



[0089] For sensor blending where / > r, it is required to find a blending matrix M e R 1 * 1 such 
that the system Sm =[A, B, MC, MD] is passive. Further, it is required to maximize the size of the 
parameter box p in which the system remains passive. Consider the parameter box defined by 



ISURF 3045 L VM 502469 



30 



{p: p, g [a. - 0S i ,u i +08,. ], / = \,2,...k} where 0 is a dilation parameter. The /th element of the 

vertex p/ j =1,2, •■•2*, of the parameter box with dilation 0 is 
given by: 



(i=l ,2 , — fc) 



(86) 



The problem is to find a blending matrix M which will maximize the dilation 0 while passifying 
the system. 



[0090] The system matrices at the vertices have the form: 



;=1 
k 



where Yy=-1 or 1. Define 

k 

A *J = -Z^M. 



i=\ 



(87) 
(88) 

(89) 
(90) 



A r ( M )P+PA(ti) PB(f4)-C T ( M )M T 
B r (/u)P-MC(/i) -(A/£>(//) + £> r (u)A/ r 



(91) 



[0091] Substituting equations 87-91 in 85,the optimal sensor blending problem can be stated 
as follows: Find matrices P e R n ' n ,{P = P T >0) and Me R r *' which will maximize 0 subject to: 



Z(p,P,M)<Q 
fory=U,...2*. 



A]jP + PA SJ PB sJ -CJjM 



B' SJ P-MC SJ -{MD SJ+ D' SJ M 



(92) 



[0092] After sensor blending, the resulting pair (MC, A) must retain observability. The 
following constraints are added to ensure observability for/? = ft: 
'A T ( M )Q 0 +Q 0 A( Al ) C T (ti)M 



<0 

[_MC{n) _i 
Qo >e 0 

where eo is the specified lower bound on Qo, a measure of observability. 



(93) 
(94) 
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[0093] Equations 92-94 define a generalized eigenvalue problem which can be solved by 
using LMI methods. If p i and p i are not known apriori, but only their nominal values (//,) are 
known, it would be necessary to guess the corresponding 5j and the maximum 0 obtained will 
depend on the choice of 8/s. 



[0094] Note that, if D =0, the off-diagonal blocks in the passivity LMI (equation 85) must be 
zero. If uncertain parameters are present only in the matrix, then robust passification is possible 
even when D =0. If, however, uncertain parameters are present in B and/or C, it will be necessary 
to have a nonzero D matrix in order to satisfy the sufficient conditions of equation 92. Thus if the 
system to be passified by sensor blending is strictly proper, an additional unknown variable D 
will have to be included in order to make the system (A,B,MC,D) robustly passive. The 
generalized eigenvalue problem is re-stated as follows for the strictly proper case. 

[0095] For the strictly proper case, the solution is to find matrices Pe B" x " ,(P =P r > 0), D 
e /F* r , and M € IF * 'which will maximize 0 subject to: 



A T ( M )P + PA(M) PB(M) ~ C T (MW T 
B T ({i)P - MC(ju) -(D + D T ) 
for/ =1,2, -2*. 



<e 



A'P + PA, PB 5j -ClM T 
B*P-MC 4 o 



(95) 



[0096] If the D matrix is added to a strictly proper system and K($)denotes the feedback 
controller for the passified system (A, B, MC, D), the effective controller (from y = Cx to w) 
becomes K eff (s) = [I + K(s)D]~ l K(s)M . If K is a constant positive definite matrix, K eff -+ D^M 
as K is increased. Thus D has the effect of reducing and limiting the effective control gain. 
Therefore, if D is a design variable, it is desirable to make it as small as possible, which can be 
accomplished by placing an upper bound on it (i.e., D +D T <d ^I). 

[0097] Note that even for the non-strictly-proper case, an additional variable D can be used 
for making (A, B, MC, MD + D ) passive. This usually increases the region of robust passivity 
but reduces the effective controller gain. 

[0098] For optimal control allocation (r>l), it is required to find a control allocation matrix N 
e R r x 'such that the system S N =[A, BN, C, DN] is passive. It is also required to maximize the 
size of the parameter box p within which the system remains passive. Define 
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ZCu,£,AO:= 



B T ( M )-C(rtX -{D{n)N 



+ N T D T (n))\ 



(96) 



[0099] The problem can be stated as follows: Find matrices S e It" * ", (S =2 T >0) and 
Ne J?""*' which will maximize 0 subject to: 

\AJZ + YAl B sj N-^C T si 
[N T Bl-C sj L — (D sj N + N T D]j) 

fory=l,2, -2* 

[0100] A lower bound constraint on a measure of controllability Q c is also needed: 



(97) 



~A{ M )Q C +Q C A T (M) B ^ )N ] <0 

_n t b t (ji) -j y 



Q C >*J 



(98) 
(99) 



where e c is the specified lower bound on Qc. 



[0101] When redundant actuators and sensors are present, it would be advantageous to 
perform both control allocation and sensor blending. If we let m < min (r,l). The problem is to 
find matrices MeJf ,x/ ,Nejr xm such that the system (A, BN, MC, MDN) is robustly passive 
(i.e., find a positive definite matrix P and matrices M y N such that 0 is maximized subject to 

A T ( M )P + PAQi) PB( M )N - C T { M )M T 1 U'P + PA. PB sj N - C^M 7 

N T B T ( M )P-MC(ju) -(MD(/u)N + N T D T (ji)M T \ < [n t BJjP~MC sJ -(MD sj N + N T D T sj M T 

(100) 

for j =1,2, — 2 , or, using the dual LMI for passivity, the problem is to find a positive 
definite matrix and £ and matrices M, N such that 0 is maximized subject to 
A{ M )Z + ^A T { M ) B(ju)N-TC T ( M ) 1 <0 p4 +42 B^N-ZC^M 7 

N t B t (ju)-MC(m)X -(MD(m)N + N t D t ( M )M t )] \N T Bl-MC sj Z — (MD sj N + N T D T sj M 7 



(101) 

for j =1,2, ■■•2*. Unfortunately, quadratic terms involving products of (M, N), (P, N) and (Z , M) 
are present in equations 100 and 101, and are therefore no longer LMIs. However, a possible 
procedure for solving this problem is to iterate between sensor blending and control allocation. 
After selecting an m < min (r,l), a full rank M matrix is arbitrarily selected and the optimal 
control allocation problem is solved, yielding an N matrix. Next, the optimal sensor blending 
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problem is solved for this fixed N, yielding a new M. The procedure is stopped when the 0s 
obtained are sufficiently close and large. 

[0102] The techniques of passification have been described. FIG. 5 shows a passivity based 
controller structure 400 for a one dimensional enclosure system 402 using the techniques 
described above. Speaker 404 is the control speaker and speaker 406 is the disturbance- 
generating speaker. Cff(s) compensator 408 denotes the transfer function of the feedforward part 
of the controller 400. C s (s) compensator 410 denotes the transfer function of the series part of 
the controller 400. Cfb(s) compensator 412 denotes the transfer function of the feedback part of 
the controller 400. N(s) and yperj(s) denote the Laplace transforms of the disturbance noise and 
the sensed noise respectively. 

[0103] FIG. 6 shows the experimentally obtained frequency-domain noise suppression 
response for the enclosure 402 for a broadband disturbance noise. FIG. 7 shows noise 
suppression in a three-dimensional reverberant structure for a three-tone disturbance noise at 
frequencies of 100 Hz, 225 Hz, and 290 Hz that was obtained via simulation. FIG. 8 shows the 
experimentally obtained noise suppression response for the structure of FIG. 7. The suppressed 
noise levels correspond to the "controller on" condition. 

[0104] Now that the modeling, passification, and controller design has been described, an 
example shall be provided. Analytically derived models (i.e., equations 37 and 42), for the 
acoustic-structure interaction dynamics in 3-D enclosures having configuration shown in FIG. 1, 
are used to obtain the finite-dimensional models for the laboratory apparatus shown in Fig. 9. 
The laboratory apparatus 500 is a rectangular wooden box with one of its surfaces being the 
aluminum plate. This surface can be set into vibrations to generate an acoustic disturbance 
inside the box using a piezoactuator attached at the center of the aluminum plate. System 
parameters for this 3-D acoustic enclosure are given in Table 1 . 



Parameter 


Value 


Parameter 


Value 


/x 


0.3048 m 


/v 


0.3810 m 


h 


0.7874 m 


a 


0.0063 m 


b 


0.2985 m 


c 


0.0063 m 


d 


0.3747 m 


h 


0.000813 m 


E 


71*10 y N/m z 




0.33 


P 


2810kg/m J 




343 m/s 


Po 


1.13 kgW 


x 0 


0.1270 m 


yo 


0.1524 m 


zo 


0.3112m 




0.1524 m 




0.1905 m 



Table 1 
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[0105] Using the parameter values of Table 1 and assumed modes approximation, a transfer 
function model of the system is obtained using equation 42. For finite dimensional 
approximation, the first five modes in each x, y, and z direction are considered, which yield a 
total of 125 modes. 

[0106] In order to validate the analytical model, an experimental frequency response was 
obtained using swept sine excitation of the piezoactuator as well as the speaker and measuring 
the acoustic pressure at the sensing microphone. The sensing microphone is located at location 
(171, 177, 285) mm. The experimental frequency response is given in Fig. 10. Using this 
experimental frequency response data, the state-space model of the system was derived using 
system identification techniques. The system identification toolbox SOCIT (developed at 
NASA LARC) was used for this purpose. The SOCIT toolbox uses ERA (Eigensystem 
Realization Algorithm) method, which is based on Markov parameters and singular value 
decomposition techniques. The ERA uses modal amplitude coherence for ranking of the most 
effective modes of the system. For the system under consideration, this identification algorithm 
gave a very satisfactory match as shown in Fig. 1 0. The identified model is a 40th order model. 
It is to be noted that identification becomes increasingly difficult as we go from SISO to MIMO 
systems. For the MIMO case, it is very hard to obtain a good match and it invariably yields very 
high order modes. 

[0107] The natural frequencies obtained from the identified system model were compared 
with the natural frequencies determined from analytical model to assess the accuracy of the 
analytical model. The comparison is given in Table 2. The symbol * in the second column 
denotes the frequencies that are missing from the analytical model but appear in identified 
model. Similarly, the symbol # in the third column denotes the frequencies that are missing in 
the identified model but appear in analytical model. The missing dynamics in the analytical 
model is attributed to the fact that the analytical model does not account for actuator, sensor, and 
other hardware dynamics and the boundary conditions used to derive analytical model are only 
approximate. In particular, the modes 2, 7, 1 1, 14, and 19 are missing in the analytical model. 
The modes 3, 9, 10, 18, 21, 22, and 24 are missing in the identified model. The modes 9 and 10 
are very close to 8 and with appropriate damping may collapse onto one. This may be the case 
for the missing modes 9 and 1 0 in the identified model. For the same reason, modes 1 8 and 24 
are missing in the identified model. The modes 3, 21 and 22 seem to be highly damped in the 
experimental data which is why the identification algorithm cannot detect. 

[0108] The frequencies listed in Table 2 do not account for boundary impedance, i.e., it 
assumes the perfectly hard boundary with zero transmissibility. For controller design, the 
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identified system model is used as the design model since it represents the system more 
accurately than the analytical model. Although the natural frequencies from analytical model 
match reasonably well with natural frequencies from experimental data, the damping on each 
mode needs to be adj usted manually since there is no systematic method to derive it 



Mode no. 


Frequency (Hz) 


Analvtical 


Kxneri m en tal 

i_/AUvl illlvlllvU 


1 


37.5 


39.2 


2 


* 


64.0 


3 


80.8 


# 


4 


106.4 


98.4 


5 


14Q K 

1 t/.O 


19S S 


6 


153.2 


1 S7 7 


7 


* 


1 89 4 


g 


91 7 8 


999 4 


9 


991 4 


TT 


10 


999 1 


TT 


11 




9^0 0 


12 


9S4 4 


9S1 1 


13 


964 8 


969 1 


14 


* 


280.5 


15 


323.4 


327.1 


16 


337.1 


351.0 


17 


382.3 


381.0 


18 


384.6 


# 


19 


* 


416.7 


20 


425.7 


423.5 


21 


435.6 


# 


22 


438.4 


# 


23 


450.1 


454.5 


24 


453.6 


# 


25 


498.0 


479.0 



Table 2 



[0109] AAs can be seen from the above description, the approach taken here is to first 
attempt designing a feedforward compensation to passify the system. If this is not acceptable 
(for example, in the case of non-minimum-phase systems, large gain matrices may be required), 
a constant-gain feedforward compensation is designed to only to render the system minimum- 
phase. The resulting system is then rendered PR by a series (pre- or post-) compensation, or if 
applicable, by sensor blending and/or control allocation. The next step is to design a PR 
controller so that the resulting closed-loop response shows the desired noise reduction. 

[0110] All references, including publications, patent applications, and patents, cited herein 
are hereby incorporated by reference to the same extent as if each reference were individually 
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and specifically indicated to be incorporated by reference and were set forth in its entirety 
herein. 

[0111] The use of the terms "a" and "an" and "the" and similar referents in the context of 
describing the invention (especially in the context of the following claims) are to be construed to 
cover both the singular and the plural, unless otherwise indicated herein or clearly contradicted 
by context. The terms "comprising," "having," "including," and "containing" are to be 
construed as open-ended terms (i.e., meaning "including, but not limited to,") unless otherwise 
noted. Recitation of ranges of values herein are merely intended to serve as a shorthand method 
of referring individually to each separate value falling within the range, unless otherwise 
indicated herein, and each separate value is incorporated into the specification as if it were 
individually recited herein. All methods described herein can be performed in any suitable order 
unless otherwise indicated herein or otherwise clearly contradicted by context. The use of any 
and all examples, or exemplary language (e.g., "such as") provided herein, is intended merely to 
better illuminate the invention and does not pose a limitation on the scope of the invention 
unless otherwise claimed. No language in the specification should be construed as indicating 
any non-claimed element as essential to the practice of the invention. 

[0112] Preferred embodiments of this invention are described herein, including the best 
mode known to the inventors for carrying out the invention. Variations of those preferred 
embodiments may become apparent to those of ordinary skill in the art upon reading the 
foregoing description. The inventors expect skilled artisans to employ such variations as 
appropriate, and the inventors intend for the invention to be practiced otherwise than as 
specifically described herein. Accordingly, this invention includes all modifications and 
equivalents of the subject matter recited in the claims appended hereto as permitted by 
applicable law. Moreover, any combination of the above-described elements in all possible 
variations thereof is encompassed by the invention unless otherwise indicated herein or 
otherwise clearly contradicted by context. 



